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We study the graviton propagator in euclidean loop quantum gravity, using the spinfoam formalism. 
We use boundary-amplitude and group-field-theory techniques, and compute one component of 
the propagator to first order, under a number of approximations, obtaining the correct spacetime 
dependence. In the large distance limit, the only term of the vertex amplitude that contributes 
is the exponential of the Regge action: the other terms, that have raised doubts on the physical 
viability of the model, are suppressed by the phase of the vacuum state, which is determined by the 
extrinsic geometry of the boundary. 



An open problem in quantum gravity is to compute particle scattering amplitudes from the 
background-independent theory and recover low-energy physics 0]. The difficulty is that general 
covariance makes conventional n-point functions ill-defined in the absence of a background. A strat- 
egy for addressing this problem has been suggested in 0; the idea is to study the boundary amplitude, 
namely the functional integral over a finite spacetime region, seen as a function of the boundary value 
of the field Q. In conventional quantum field theory, this boundary amplitude is well-defined (see 
0]) and codes the physical information of the theory; so does in quantum gravity, but in a fully 
background-independent manner |^. A generally covariant definition of n-point functions can then 
be based on the idea that the distance between physical points -arguments of the ri-point function- is 
determined by the state of the gravitational field on the boundary of the spacetime region considered. 
In this paper we implement this strategy. 

We compute a first order term of the (connected) two-point function, starting from full non- 
perturbative quantum general relativity. To this order, the computation produces a 4d version of 
the "nutshell" 3d model studied in [^. Using a natural gaussian form of the vacuum state, peaked on 
the intrinsic as well as the extrinsic geometry of the boundary, we derive an expression for a compo- 
nent of the graviton propagator. At low energy, this agrees with the conventional graviton propagator. 
Other components and second order terms will be presented elsewhere 0. Our main motivation is 
to show that a technique for computing particle scattering amplitudes in background-independent 
theories can be developed. (The viability of the notion of particle in a finite region is discussed in 0. 
For the general relativistic formulation of quantum mechanics underlying this calculation, see |^ . On 
the relation between graviton propagator and 3-geometries transition amplitudes in the conventional 
perturbative expansion, see [lO|-) 

We use some standard loop quantum gravity (LQG) results as well as a specific spinfoam 

model 0,0: the GFT/B theory, in the terminology of 0, defined using group-field-theory methods. 
The results do not change using the GFT/C model. These are tentative background-independent 
quantizations of euclidean general relativity. The first was introduced i n 1131 and is favored by some 
recent arguments 0, 0| . The second was introduced in |0 (see also |l7l|') and has good finiteness 
properties 

The physical correctness of these theories has been questioned because in the low-energy limit their 
interaction vertex (lOj symbol, or Barrett-Crane vertex amplitude) has been shown to include -beside 
the "good" term approximating the exponential of the Einstein-Hilbert action also two "bad" 
terms: an exponential with opposite sign, giving the cosine of Regge action [ijj (analogous to the 
cosine in the Ponzano-Regge model) and a dominant term that depends on the existence of degenerate 
four-simplices |2^0|. We show here that only the "good" term contributes to the propagator. The 
others are suppressed by the rapidly oscillating phase in the vacuum state that peaks the state on its 
correct extrinsic geometry. Thus, the physical state selects the "forward" propagating component 
of the transition amplitude. This phenomenon was anticipated in |23j . 
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Following the 2-point function of quantum general relativity can be extracted from a nonper- 
turbative formalism using the expression 

W^'^^ix, y; q) WW] {s'\h^''{xW{y)\s) (1) 

ss' 

Here q is the intrinsic and extrinsic (a coherent state depends on both classical position and momen- 
tum) geometry of a closed 3d surface, interpreted as the physical boundary of a 4d spacetime region, 
and X and y are points on this (metric) surface. W[s] is the boundary amplitude, that codes the 
quantum gravity dynamics, ^'g[s] is a quantum state of the geometry peaked on q, the sum is over a 
basis s of states of the 3-geometry, namely s-knot states, and h°-^{x) is a linearized gravitational field 
operator acting on the state vE'q(s)- We refer to ^] for a motivation and a discussion of this expression 
and to for notation and an introduction to the formalism and its ideas. TV is the normalization 
factor defined by 1 = A/'X^^ W^[s]^'g(s). Eq is well-defined if we have a dynamical model giving 
W[s], a "vacuum" state vE'q[s] and a form for the operator h°''^{x). 

We choose the boundary functional W\s\ defined by the group field theory GFT/B. We recall the 
definition of W[s], referring to ^ and l^j for motivations and details. The theory is defined for a 
field (p '■ (*S'0(4))'* ^ i? by an action of the form 

S[cf>] = S^iM + ^SinM. (2) 

The field (/) can be expanded in modes (l>a^[['',a^i (Eq. (9.71) of |^). Indices j„, n = 1, 4 label simple 
50(4) irreducible representations. Irreducible representations of SO{4) are labelled by a pair of spins 
j_), corresponding to the split of so(4) = su{2) x su{2) into its self-dual and antiself-dual rotations; 
simple representations are those where j-^- = j- = j, and are therefore labelled by a single spin j. The 
index q;„ labels the components of vectors in the representation j„. The index i labels an orthonormal 
basis of intertwiners (invariant vectors) on the tensor product of the four representations j„. We 
choose a basis in which one of the basis elements is the Barrett-Crane intertwiner igc (Eq. (9.99) of 
1^). S'0(4)-invariant observables of the theory are computed as the expectation values 

W[s] = I D0 fM e-/^'-*/*', (3) 

where fs{4>) is the function of the field determined by a spinnetwork s = (T^ji^in). Recall that 
a spinnetwork is a graph F formed by nodes n connected by links I, colored with representations 
ji associated to the links and intertwiners in associated to the nodes. We indicate by Inm a link 
connecting the nodes n and m and by j„m = jmn the corresponding color. The spinnetwork function 
is defined by 

fsW-EIKlZ^.-^ (4) 



here n runs over the nodes and, for each n, the index m runs over the four nodes that bound the 
four links Inm joining at n. Each index anm = ctmn appears exactly twice in the sum, and is thus 
contracted. The perturbative expansion of l(2Jl in A leads to a sum over Feynman diagrams. The 
Feynman rules are given by the propagator 

= kr' E n ^^^A,., 5^r.<.U.V (5) 
ir(n) n 

where 7r(n) are the permutations of the four numbers n — 1, 2, 3, 4; and the vertex amplitude 

KT^" = ( n *»"«c) ( n ^".^"..0 l^ijnrn), (6) 

n n<7n 
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where the index n = 1,...,5 labels the five legs of the five-valent vertex; while the index m ^ n 
labels the four indices on each leg. B{jnm) is the lOj symbol, or Barrett-Crane vertex amplitude 
[l^. whose large j expansion is discussed below. The sums over permutations in the propagator 
give rises to a number of terms. Each of these can be interpreted as a spinfoam cr, by identifying 
closed sequences of contracted deltas as faces. Hence the amplitude lj2Jl can be written as a sum of 
amplitudes of spinfoams bounded by a given spinnetwork W[s\ — X](5o-=s ^["'l' expression that 
is naturally interpreted (and can also be derived) as a sum over discretized 4-geometries bounded 
by a given discretized 3-geometry, namely as a definition of the Misner-Hawking sum-over-geometries 
formulation of quantum gravity, where the amplitude of a 3-geometry ^5 is given by the Feynman 
integral over 4-geometries W^g] = Jgg^3gDg e^'^'^l. The amplitude A[a] of a spinfoam a with n 
vertices is proportional to A". Below we consider the first order of the expansion in A. Since each 
face / gives then a contribution dimj/ (by contracting the deltas), in the large j limit the dominant 
term is the one with the largest number of faces. To first order in A, we have only one possibility: one 
vertex and a five-valent boundary spinnetwork. The dominant contribution for large j is obtained by 
choosing the spinfoam a and the spinnetwork s = da defined by the dual of a four-simplex and its 
boundary. The spinnetwork s has five nodes, connected by 10 links as follows 

il 



(7) 




24 j34 *3 

The boundary function /s((^) determined by this spinnetwork is fs{4>) = X]q„„ nri=i 5 '^iTrT*"- This 
is a an observable in the group field theory. The dominant term of its expectation value jSj is 

= ^ (n(*"l*Bc) j [ n dim(i,„„) j BUnra). (8) 

The second ingredient we need is a boundary state ^'(j[s]. To identify it, we need the geometrical 
interpretation of the boundary spinnetwork s. To this aim, we use the fact that the spinfoam model 
can be obtained from a discretization of general relativity on a triangulated spacetime. Introduce 4d 
coordinates and represent the gravitational field by means of the one- form tetrad field (x) = 
ej^{x)dx^ (related to Einstein's metric by g^^{x) = e^(a;)e/^i(a;)). Assuming that the triangulation 
is fine enough for this field to be approximately constant on a tetrahedron, with constant value e^, 
associate the 4d vector = e^Ax^ to the segment s of the triangulation, where Ax^ is the coordinate 
difference between the initial and final extremes of s. To each triangle t of the triangulation, associate 
the bivector (the object with two antisymmetric indices) 

Bl' ^eiei-eiel, (9) 

where s and s' are two sides of the triangle. (As far as orientation is kept consistent, the choice of the 
sides does not matter). Bl'^ is a discretization of the Plebanski two-form B^'^ — l\e'^ . The quantum 
theory is then formally obtained by choosing the quantities B\'^ as basic variables, and identifying 
them with 5*0(4) generators J/"^ associated to each triangle of the triangulation, or, equivalently, to 
each face of the corresponding dual spinfoam. The form © implies that tijKhBl^ B^^ = any time 
t = t' or t and t' share an edge. Accordingly, the pseudo-scalar Casimir C = eijKLJt"^ Jf"^ = is 
required to vanish. This determines the restriction to the simple representations, which are the ones for 
which C — 0. The scalar Casimir C — ^J{'' Jtij — ^Bl'-' Btij, on the other hand, is easily recognized, 
using again lO, as the square of the area At of the triangle t. For simple representations, the value 



4 



of C is j{j + 1). The quantization of the physical area, with eigenvalues proportional to + 1) is 

a key result of LQG reappearing here in the context of the spinfoam models. The LQG result 
assures us that we can interpret it as a physical quantization and not an artifact of the discretization 
and fixes the proportionality constant: Aj — ^ttTiG^J j{j + 1). The geometrical interpretation of the 
intertwiners can be obtained in a similar way and will be discusses elsewhere. 

Let now q be the intrinsic and extrinsic geometry of the boundary Yi„ of a 4d ( metric ) sphere. We 
want to construct the state ^g[s]. (On the vacuum states in LQG, see |2^l25lEa l2^Eq.) Below we 
shall only need the value of ^ q[s\ for the spinnetworks defined on a graphs F, dual to 3d triangulations 
A. We identify each such A with a fixed triangulation of Sg. We assume here for simplicity that 

q{s\ is non vanishing only for s whose intertwiners are ibc intertwiners. This choice preserves the 
symmetry of the sphere. (There is another simple symmetric choice, used in the first version of this 
draft, which is to choose '^q[s\ independent from the intertwiners. This possibility will be explored 
elsewhere.) The areas Anm of the triangles tnm of A determine background values jnm of the spins, 
via = (SnhG)^ jnmijnX + 1). We want a state ^q[s] = ^q{T,jnm) peaked on these background 
values. The simplest possibility is to choose a gaussian peaked on these values. However, this leaves 
the possibility open of having a phase 

^q[s] = exp I -^a(^„rn){pq) {jnm ^ j^ni){jpq - j^f) + « ^nljnrn \ (10) 

where {n,m) runs on links of s, a(^nm){pq) is a numerical matrix that we will fix later on and the sum 
Sn<mp<(j i^ understood. The phase factor in this state is important. As we know from elementary 
quantum mechanics, it determines where the state is peaked in the variables conjugate to the spins 
jmn- Recall the form of the Regge action S'Rcggc = J2n<m ^nmijmn)jnm, where $„m(jm„) are dihedral 
angles at the triangles, which are function of the areas themselves and that dSjicggc/djnm = ^nm- It 
is then easy to see that these dihedral angles are the variables conjugate to the spins. Notice also that 
they code the extrinsic geometry of the boundary surface, and in GR the extrinsic curvature is the 
variable conjugate to the 3-metric. The value of ^'^[s] on the five-valent spinnetwork |(7|) can be deter- 
mined by triangulating Eg with the 3d triangulation formed by the boundary of a regular four-simplex 

of side L. The area of the triangles is A]^ — \/3i^/4. Then jmn = where %TThGyJjL{3L + 1) = Al- 
In the large L limit we take = STrhGA^. The dihedral angles ^nm = $ of a regular tetrahedron 
are given by cos(<I>) = —1/4. The 1/jl dependence of a (absent in the first version of this paper) 
ensures that the relative uncertainties of areas and angles -that is: intrinsic and extrinsic geometry 
of the boundary- become small in the large distance limit. Its need has been pointed out by Simone 
Speziale in the 3d context 29] and by John Baez in the 4d case, following numerical investigation by 
Dan Christensen and Greg Egan, that have shown that in the absence of this dependence the width of 
the gaussian is not sufficient for the approximation taken above to hold |30| | . To respect the symmetry 
of the symmetry of the sphere, the covariancc matrix a(^nm){pq) of the gaussian can depend only on 
three numbers 

(11) 

where 5(„m)(pg) = 1 if (nm) = (pq), ai^nm){pq) = 1 if just two indices are the same, and b(^nm){pq) = 1 
if all four indices are different, and in all other cases these quantities vanish. (We will use this 
notation, namely a(i2)(i3) = cti, Q;(i2)(i2) = ct2, Q!(i2)(34) = Q^a repeatedly.) The state (|10|l is thus 
completely determined up to the three numbers ai, a2, 0:3. This amounts to select a vacuum state 
which is a coherent state peaked both on the background values of the spins (the intrinsic geometry 
of the boundary surface), and on the background values of the angles (the intrinsic geometry of the 
boundary surface). See ja| for a similar construction in 3d. 

The third ingredient we need is the graviton field operator. This is the fluctuation of the metric 
operator over the flat metric h°-^{x) = g"'''{x) — S"-^. It is more convenient to consider here the 
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fluctuation of the densitized metric operator h°'''{x) = {det g)g''''{x) - 5"^ = E''\x)E^'{x) - (5"^ In 
the hnear theory, the propagators of the two agree because of the trace-free condition. To determine 
its action, we can equally use the geometrical interpretation discussed above, or LQG. We identify the 
point X with one of the nodes n of the boundary spinnetwork s. Equivalently, with (the center of) one 
of the tetrahedra of the triangulation. Four links emerge from this node. Say these are e/, / = 1, 2, 3, 4. 
They are dual to the triangular faces of the corresponding tetrahedron. Let be the oriented normal 
to this face, defined as the vector product of two sides. Then = i?"'(x)n^ can be identified 

with the action of the SU (2) generator J' on the link e/. We have then immediately that the diagonal 
terms define diagonal operators 0, |0| 

E'\n)El{n)\s) = {SirhGf + l)\s), (12) 

where ji is the spin of the link /. Non-diagonal terms, that act on the intertwiners, are discussed 
elsewhere. 

We have now all the elements needed to compute the expression Inserting ijSJ, (|10|l and H12(l 
into Q we obtain a well-defined expression for the propagator. We choose the points x and y to 
be two distinct nodes of the boundary spinnetwork. Equivalently, these can be thought as points 
located, say, in the centers of the corresponding dual tetrahedra: in the theory, of course, position is 
not determined with precision lower that the individual "atoms of spaces" described by the individual 
tetrahedra. We consider the ten by ten matrix W{L) formed by the "diagonal" components of the 
propagator 

W^(X)(,,)(,i) = W'^^ix, y- q) naUb n^n^, (13) 

where n = riij and n = and is the normal to the triangle Uj. We also write W{L) = W{L)l\n\'^. 
By symmetry 

VF(L)(y)(fe;) = Wi{L) a(ij)(^ki) - W2{L) %)(fe;) + Wz{L) 6(y)(fe;). (14) 

Before computing these quantity in the background independent theory, let us compute it in conven- 
tional linearized quantum general relativity. In a fiat background metric, two points in the center of ad- 
jacent tetrahedra, in a surface with the boundary geometry chosen, are at a distance jx" — x'l\q = L/4. 
If the four indices i, j, fc, I are all distinct, it is easy to see that n and n are orthogonal; then the prop- 
agator is easily computed to be [lo| 



.SvrTiG 1 .32?iG 

„d\2 



TT^lincarizcd / 7- \ .v^fi(</v^ _^j^,u^ /i rr\ 

iL)-^ U«_,^|2-*^ (15) 



On the other hand, the components a-^id ^(i^)(ik)'^'^ ^-''^ vacuum expectation values at fixed 

"time" : the first is the fiucuation of the area square of a triangle, and the second is the vacuum corre- 
lation between the fluctuations of the area squares of two adjacent triangles in the same tetrahedron. 
These are also proportional to L~^. We can therefore write 

^linearized = W (16) 

where W is a numerical matrix, with the same symmetry structure as in (I14|l . In particular, yV(i2)(34) = 
Wa — i, while Wi and W2 are real numbers of the order of unity, easily obtained from the linear theory. 

We now compute the matrix W{L) in the full theory. Since this is a diagonal term in the propagator, 
we can use (fT^ and (Q) reads 

W{L)(,^)(ki) = E^W ((8^nG)2j-,(j,, + 1) - \n\^mi:hGfju{jM + 1) - \n\'')^M- (17) 
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The terms |np come from the background S""^ and are equal to the square of the area of the face, 
namely to {SnhGjL)^ , for large j. Inserting ((SJ and (ITUIl we have 

W?(L)(,,)(,;) = J^h(fiipl ^ I Yl dim(j„„0 J + 1) - jI) {jkiijki + 1) - jI) 

jnm \n<m / 

B{jnm) exp \ -^a(nm){pq){inm ~ jL){jpq - j l) + i'^^jnm \ (18) 
I ^-^^ n,m J 

Since the vacuum exponential peaks the sum around jl, which is large, we can discard the +1. We 
expand the summand in the fluctuations 5jij = (jy — jl), and keep only the lowest term, assuming 
that the gaussian suppress the higher terms. We redefine N absorbing all the factors that cancel 
with terms in (the expansion to this order of) the partition function. We assume that the dimj terms 
vary slowly over the range where the gaussian is peaked, and can be considered constant. We change 
summation variable from the spins to the fluctuation of the spins 

W?(L)(,,)(fe,) = AN{SnhGfjl J2 Sj.j SjM B{jl + Sj^m) e-W"(-)(-)^-^-*^-+^*E_J"'". (19) 

The rapidly oscillating term expi^^^^^j^ j„,„ tends to suppress the sum. To evaluate it, we need the 
explicit form of B{jL + 5 jnm) in the large j regime. This is of the form [20ll2l| 



B{3nm) = Y.P- [e'^'^--^"^+'=-^ + e-'^^--(")+'=-^] + D{jnm) (20) 



Here r labels the distinct 4-simplices having areas A,„„ = jnm {jnm + 1) and is a slowly varying 
factor. Since the sum H19|l is peaked around jnm = jL, let us expand the 10 j symbol around this 
point. To second order around jnm = jL, the Regge action reads 

*S'Rcggc(jn7n) ^ ^ ^ Jnm ^ (^nin){pq)^jmn^jpq, (21) 
nm 

where, introducing the "discrete derivative" = /(j + 1/2) — /(j), we have defined 



d^mnijrs) 



G(mn)(pq) q-^^ 



(22) 

=3L 



Thus, around jnm — jh, we have 



B {jnm ) 



+ D{jnm) 

(23) 

where is the regular four simplex (for which = 1), which is only non-degenerate four-simplex 
with these areas |20|. The key observation is now the fact that the rapidly oscillating exp{i$ X^nm jnm} 
term in the second term of this expression cancels with the rapidly oscillating term in p9|l . Therefore 
the second term of the last expression contributes in a non- negligible way to the sum H19() . The first 
term is suppressed (by the rapidly oscillating factor exp 2i$ j„m) and it is reasonable to expect 
that so is the degenerate term D{jnm) because this corresponds to 4-simplices with different angles, 
and should be dominated by different frequencies. Therefore 1)19(1 becomes, keeping only the first term 

iy(L)(,,)(feO - 4J^j fP{rR)e'^ ^ 6j,, Sjm e*G(„™,(,,)5,„„.5,,,g-5fe"(.-.)(P,)^J-"."*>,^ (24) 
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where we have also used W{L) ~ \n\^W{L) = {8TTfiGjL)^W{L). Approximating the sum by a gaussian 
integral gives 

W{L) ^Aj^^jZ' a -zCy' (25) 

We only need to evaluate the derivatives (|22ll of the angles with respect to the spins. 4d geometry 
gives |3ll | 

"dA ~ \/hL? (^2*^*-"™'"'^'-' ~ ^'^("™)(p<?) ^ 46(„m)(p^)^ = — /C(„m)(pg). (26) 

The ten by ten matrix /C has purely numerical entries. From the relation between areas and spins, we 
have 

G(nm)ipq) = S,T:hGL~^ IC(^nm){pq) = ^{4:] l}^^ )^(nm)(pq) ■ (27) 

Notice the Jl factor that combines with the one in front of a in H25I) to give the crucial overall 1 /Jl 
dependence of the propagator. Using the rcation between spins and areas, 125|l reads 

This is precisely the value p5ll6|l of the propagator computed from the linearized theory, with the cor- 
rect l/\x — spacetime dependence. The three numerical coefficients of the matrix a are completely 
determined by a = 47rVV3 W'^ - i%/3/4 /C. 

The propagator we have computed approximates the perturbative one. This is correct for large 
|x — yl, where the approximations used in evaluating the sum IjlSI) hold. When \x — y\ approaches the 
Planck length, quantum gravitational corrections appear. In this regime, H18() can be easily computed 
numerically. Given the discreteness of the representation, the propagator is likely to be cut-off at the 
Planck scale. Therefore we obtain the result that, at least to this order, the (component we have 
computed of the) propagator is the standard free-theory one at large distances, with corrections at 
the Planck scale, wilch appear to make it finite. 

Our main results are: (i) the "bad" terms of the Barrett-Crane vertex amplitude are cancelled 
when the correct vacuum state, appropriately peaked on the extrinsic geometry, is used; (ii) the 
correct propagator emerges to this order; and (iii) --for us the most interesting result- it appears to 
be possible to compute n-point functions from background-independent quantum field theories. 

I am strongly indebted to Leonardo Modesto, Simone Speziale, John Baez and Dan Christensen for 
crucial criticisms, inputs and suggestions. 
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